Spectral properties of the dimerized and frustrated S = 1/2 chain 
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Spectral densities are calculated for the dimerized and frustrated S = \ chain using the method 
of continuous unitary transformations (CUTs). The transformation to an effective triplon model is 
realized in a perturbative fashion up to high orders about the limit of isolated dimers. An efficient 
description in terms of triplons (elementary triplets) is possible: a detailed analysis of the spectral 
densities is provided for strong and intermediate dimerization including the influence of frustration. 
Precise predictions are made for inelastic neutron scattering experiments probing the S = 1 sector 
and for optical experiments (Raman scattering, infrared absorption) probing the 5 = sector. 
Bound states and resonances influence the important continua strongly. The comparison with the 
field theoretic results reveals that the sine-Gordon model describes the low-energy features for strong 
to intermediate dimerization only at critical frustration. 
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I. INTRODUCTION 

One-dimensional quantum antiferromagnets display 
very interesting and fascinating physical properties. The 
generic system to be studied in this context is the dimer- 
ized and frustrated Heisenberg chain which comprises 
many physical phenomena of interest. There are gap- 
less and gapped phases^, fractional excitationsS*^ and 
confinement^. 

The quantitative calculation of spectral densities is a 
very important issue in the field of strongly correlated 
electron systems. The interplay between kinetics, inter- 
action and matrix elements leads to characteristic struc- 
tures in the spectra. In recent years there has been signif- 
icant progress in calculating spectral densities in the field 
of quasi one-dimensional quantum spin systems, e.g., the 
uniform Heisenberg chaini, spin ladders systems^ and 
strongly dimerized spin chainsL^. 

Beside the fascinating theoretical aspects, the de- 
termination of spectral densities is of direct impor- 
tance for experimental measurements. The theo- 
retically predicted spectra are relevant for inelas- 
tic neutron scattering experiments and optical ex- 
periments like Raman spectroscopy and infrared ab- 
sorption. There is a large number of quasi one- 
dimcnsional compounds which can be successfully de- 
scribed by the dimerized and frustrated Heisenberg 
model, e.g., spin-Peierls compounds like CuGeO j^d°d^d^ 
and a'-NaV203i^, (VO)2P207 — or organic compounds 

,)2 . 2.5D20)iS*ii, 
compounds like 



like Cu2(C5Hi2N2)2Cl4^^ and (Cu(N03)2 ■ 2.5D20)iS*ii, 



112iN2;2'-^i4 

but also undimerized cuprate chain 
KCuFgia, Sr2Cu05i9 and SrCnO-^. 

We will describe the dimerized and frustrated spin 
chain in terms of elementary triplets (triplons)^^ which 
carry total spin one. The commonly accepted elementary 
excitations for the undimerized system are fractional ex- 
citations, so-called spinous carrying S = 1/2^'^. Recently 
we have shown that a description in terms of triplons is 
also possible for the isotropic Heisenberg chain^^. So 



there is no necessity to use fractional excitations in one- 
dimensional systems. Remarkably, even more spectral 
weight is captured by the states of two triplons than 
with two spinous^. This result was recently confirmed 
by Hamer et al&. 

In this work, we take a closer look at the generic fea- 
tures of the spectral properties at finite dimerization. In 
this regime a description in terms of triplons is surely cor- 
rect due to the confinement of the spinous^. We will learn 
more about the triplon picture and try to connect the 
triplon-based findings at finite dimerization with results 
obtained at zero dimerization using the spinon picture. 

The paper is organized as follows. Sect. II gives an 
introduction in the model we consider in this work. 
Sect. HI presents the method we use and introduces the 
basic quantities under consideration. In Sect. IV results 
for the dynamical structure factor are shown. We present 
results for the one-triplon contribution, the two-triplon 
contribution and compare with field theoretical results. 
Sect. V shows spectral densities with total spin zero. We 
give results for the two-triplon contribution for nearest 
neighbor (NN) and next-nearest neighbor (NNN) cou- 
pling. We also provide results for Raman spectroscopy 
and optical absorption. Sect. VI summarizes this paper 
and Sect. VII comprises the conclusions. 



II. MODEL 

The Hamiltonian for the dimerized and frustrated S = 
1/2 spin chain reads 

i/ = Jo ^ [(1 + <5(-l)^) S,S,+i + aoS,S,+2] , (1) 



where 5 parameterizes the dimerization and ao the rela- 
tive frustration between next-nearest neighbor spins. In 
order to apply a perturbative treatment we transform 
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FIG. 1: Phase diagram of the dimerized and frustrated 
Heisenberg chain depending on frustration ao and dimeriza- 
tion S. The system is always in a gapped regime except in 
the interval ao £ [0, ao,c] at zero dimerization. The dashed 
line marks the Shastry-Sutherland line where the ground-state 
is known exactly. Solid squares correspond to the parame- 
ters (1 - S)/{1 + 5) ^ X € {0.3; 0.6} and ao/{l - S) = a € 
{0.0; 0.25; 0.5} used in this work. 



Eq. (QJ into 

H/J = ^ [821821+1 + AS2iS2i-i + Aa8i8i+2] , (2) 

i 

where J = Jq{1 + S), X = {1 - S)/{1 + 6) and a = 
ao/{l-6). 

The dimerized and frustrated spin chain exhibits very 
interesting intrinsic physics. The phase diagram of the 
model is shown in Fig.^ At 6 — there are two regimes. 

(i) For ao < ao,c the ground state is in the same uni- 
versality class as the uniform Heisenberg chain. The 
excitations are massless and the standard description 
is in terms of unconfined spinous carrying total spin 
S = 1/2^ '^. In a recent work we have shown that also 
a description in terms of triplons is appropriate for the 
isotropic Heisenberg chain^S,. 

(ii) At ao = ao,c there is a transition into a spon- 
taneously dimerized phassiSiSiSiiS^. The ground-state 
is two-fold degenerate and the excitations are massive 
spinous for ao > ao,c- At the Majumdar- Ghosh point 
(ao — 0.5) the ground-state is known exactljs2Si2L2^. The 
validity of a triplonic description in this gapped phase is 
an open issue. On the other hand, Haldane^ has shown 
that for any finite dimerization S the spinous become 
confined and a description in terms of triplons carry- 
ing total spin 5' = 1 is appropriate. The spectrum is 
always gappedi*22i2& and the excitations can be viewed 
as bound states of two spinon a^^'^^i"^^'^"^ . The interac- 
tion between the spinons is not exhausted by the con- 
finement and there remains a triplon-triplon interaction 
which can lead to two-triplon bound states with total 



spin S = and S* = 1 lying below the multi-triplon 
cont inuurr^SiSiSiSSiSLSi^SiiMi 

Spectral properties are particularly difficult to address. 
So far, results can be obtained either by numerical ap- 
proaches like exact diagonalization or quantum Monte 
Carlo or by studies of effective continuum models. In 
particular, the case of finite frustration and the 5 = 
sector relevant for optical experiments has not yet been 
investigated thoroughly. 

In the following we expand about the limit of isolated 
dimers on the strong bonds, i.e A = 0. We present results 
for the two-triplon contribution to the spectral density for 
strong (A = 0.3) and intermediate (A = 0.6) dimerization 
and for three representative values of the frustration (a = 
0,a ~ 0.25, a — 0.5). These parameters are marked in 
Fig. n as filled squares. 



III. METHOD 

A continuous unitary transformation-^ is used to map 
the Hamiltonian H to an effective Hamiltonian H^s 
which conserves the number of triplons on the strong 



bonds, i.e. [Ho,Hcff] — where Hq := H\ 



43.46 



The 



ground state of Hcs is the triplon vacuum. Using an 
infinitesimal antihcrmitean generator 77 we have 



(3) 



where H{0) ~ H, H{oo) — Hcs and I is an auxiliary 
variable. The optimized choice for 77 reads 



(4) 



where the matrix elements rjij and Hi j are given in an 
eigen-basis of Hq counting the number of triplons on the 
strong bonds. The choice Q) retains only triplon conserv- 
ing processes and it eliminates all parts of H changing 
the number of triplons^^. It must be stressed that the 
approach represents a renormalization in the sense that 
matrix elements between states with very different ener- 
gies are transformed strongest. The matrix elements be- 
tween energetically similar states are transformed only at 
a later stage of the transformation. In order to determine 
spectral weights /n and spectral densities the observable 
R is mapped to an effective observable by the same 
unitary transformation^^ ^ as the Hamiltonian. 

The continuous unitary transformation cannot be car- 
ried out without truncation. We use a perturbative 
method in A. The effective hamiltonian is calculated 
up to order 10 and the effective observable is calcu- 
lated up to order 7 in the two-triplon sector. The plain 
series of the important quantities will be made available 
on our home pagea^. The numerical effort to obtain the 
result is discussed in detail for the spin-ladder systems 
and will not be discussed in this work^. 

The following extrapolation technique is employed 
which we have introduced recently for spin-ladder 
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systemsaiaa. After fixing a to the desired value the plain 
series in A is converted into a series in 1 — A invoking the 
one-triplon gap. The one-triplon gap is the natural inter- 
nal energy scale of the problem. Since in this work we are 
interested only in strong and intermediate dimerization, 
no further extrapolation techniques like standard Pade 
extrapolants are used. There is no uncertainty in the ob- 
tained spectral densities for strong dimerization A = 0.3. 
The uncertainty is about 2% for the worst case at inter- 
mediate dimerization (A = 0.6). In order to investigate 
the cases of weak or vanishing dimerization it would be 
important to use further extrapolation tools and to treat 
processes with longer or infinite range explicitly. 

Spectral densities can be calculated at T = by the 
retarded resolvent 



= --lm( 



1 



{H - Eo) 



iO- 



-R 



(5) 



Due to the conservation of triplons after the unitary 
transformation the spectral density can be split into ad- 
ditive parts In{oj), the n-triplon contribution to the spec- 
tral density, 



I{U0) = Y.In{u 



(6) 



Im( 

TT 



Rl 



iO- 



-R, 



cff ,n 



The total intensity can be calculated using Dirac's iden- 
tity for Eq. ((SJl 



dLuI{uj) = (0|i?t,<=ffi^cff|Q^ _ 



(7) 



0+ 



The total intensity /tot is a sum over the spectral weight 
of all triplon sectors, /tot = X^n^"- ^he spectral weight 
/„ in the n-triplon sector is calculated by 



/„ = (0|/?t^'=«/?f |0) 



(8) 



where Rn denotes all n-triplon excitation processes of 
the observable R. Using the sum rule /tot = (0|/?^|0) — 
(0|/?|0)^ we can check the reliability of the perturba- 
tive results. For later use we define the relative spectral 

weights /n^icl = /n//tot with J^n^n/hot = 1- 



the two-triplon channel on passing from strong to zero 
dimerization^^. In an analogous procedure we have an- 
alyzed the spectral weight distribution on the disorder 
line at a = 0.5^. We find indications that again the 
two-triplon contribution is the dominant one, even for 
vanishing dimerization. But due to the complexity of 
the frustrated system no unambiguous extrapolations are 
possible"*^. Thus we do not have a final answer for the 
massive frustrated phase. 

In the following we will show results for the one-triplon 
and the two-triplon contribution to the spectral density. 
The one-triplon contribution contains most of the spec- 
tral weight at strong and intermediate dimerization. In 
the limit of small dimerization it is highly reduced and 
becomes unimportant for zero dimerization. We expect 
that the two-triplon contribution is the only relevant con- 
tribution in the limit of zero dimerization. 

A discussion of the one-triplon contribution to the dy- 
namic susceptibility of a dimerized chain without frustra- 
tion can be found in a work by Miiller and MikeskaiS. Re- 
cently Zheng et al. published results for the one- and two- 
triplon contribution of a strongly dimerized spin chain 
without frustration^. Our results at zero frustration 
agree with the findings of Zheng et al.. 

Here we want to extract the generic features of the two- 
triplon contribution for various dimerization and frustra- 
tion in order to gain insight in the evolution of this contri- 
bution in the limit of vanishing dimerization. Therefore 
it is interesting to compare our results with results for 
the dynamical structure factor for zero dimerization. 

At a = an exact calculation of the two-spinon con- 
tribution to the dynamical structure factor using Bethe- 
ansatz is possible'*. The two-spinon contribution ex- 
hausts 72.89% of the total spectral weight and it displays 
a singular divergent behavior at the lower edge of the two- 
spinon continuum. At finite frustration only numerical 
results using exact diagonalization at finite temperatures 
including frustration are available^*. In addition, there 
are also results using abelian bosonization extracting the 
universal features of the dynamical structure factor at 
low energies for small dimerization?'*. In the following 
we will identify the major features of these studies in our 
triplonic description at finite dimerization. 



A. One- Triplon Contribution 



IV. S=l EXCITATIONS 

This part of the paper contains results for the dynam- 
ical structure factor of the dimerized and frustrated spin 
chain. At A = the system consists of isolated dimers 
and therefore the total spectral weight is in the one- 
triplon channel /i. Turning on A will reduce the spec- 
tral weight in the one-triplon channel and the spectral 
weight will also be distributed in the multi-triplon chan- 
nels. In the unfrustrated chain the spectral weight is 
shifted almost totally from the one-triplon channel in 



The physical observable R^ * for total S' = 1 excita- 
tions is locally 

Aoc ^ = Si (9) 

where i denotes a site of the chain. In the following we 
will name the left and the right site on a strong bond 
by L and R, cf. Fig. 12 We decompose the action of the 
full observable /?'^^*|itrp in the one-triplon channel on 
the ground-state for fixed one-triplon momentum k by 
writing 

/?^=^|it,p|0) = ^ft7p(fc)l^>, (10) 
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FIG. 2: Sketch of the local observable for S=l excitations. 
Double lines denote strong bonds and single lines weak bonds. 
The observable couples to the left (L) and to the right (R) 
spin of a strong bond. 



where the amplitudes Aftrp are given by 



aS=i 

^Itrp 



V2i^af sin (k{2l + ^ 



(11) 



The sum runs over all strong bonds I. The coefficient af 
is the ampHtude for the creation of one triplon at site / by 
SI- The amplitudes of Sf^ do not need to be calculated 



due to the symmetry 



The basic unit length 



a/2 is the distance between two neighboring sites. In the 
figures momenta between and 1 are given in units of 
7r/(a/2). So the comparison between our results and the 
conventional notation for undimerized chains using the 
distance a' between two neighboring spins as unit length 
(a' = a/2) is simplified. 

In Fig.|21the results for the one-triplon dispersion uj{k) 
(left panels) and the one-triplon spectral weight Ii{k) 
(right panels) are shown. We present results for A = 
{0.3; 0.4; 0.5; 0.6} and a = (Fig.OJi), a = 0.25 (Fig.Ob) 
and a = 0.5 (Fig.Et). 

Let us discuss first the case of vanishing frustration 
(a = 0). At A = the system consists of isolated dimers 
and the one-triplon dispersion is fiat. Turning on A the 
triplon starts to hop from dimer to dimer and it acquires 
a finite dispersion (Fig. left panel). The dispersion 
has minima at fc = and fc = tt (in units of l/(a/2)) 
which represent the one-triplon gap A. In the limit of 
zero dimerization the one-triplon gap closes and it is to be 
expected that the one-triplon dispersion equals the well 
known Cloizeaux-Pearson^i dispersion relation u;cp(fc) = 
7r/2|sin(fc)|. 

The one-triplon spectral weight Ii{k) is shown in the 
right panel of Fig.l^i. The leading term of Ii{k) is pro- 
portional to sin^(fc/2). It is called the dimer structure 
factor— (note again that we measure the momentum k in 
units of l/(a/2)). The one-triplon spectral weight is con- 
centrated at fc = TT. At finite dimerization the reduction 
of Ii due to the inter-dimer exchange is achieved by mo- 
menta k < 0.97r. The spectral weight increases in a small 



interval about k = tt— . In the limit of smaller dimeriza- 
tion the one-triplon spectral weight becomes more con- 
centrated about k — n. The total weight, integrated over 
momentum and frequency, vanishes for vanishing dimer- 
ization. 

In Fig. |3h> the corresponding results for a = 0.25 are 
shown. The one-triplon dispersion is similar in shape to 
the case of zero frustration. Due to the finite frustration 
the excitations become more local and the triplon is less 
dispersive. The gap values are slightly larger and the 
maximum values of the one-triplon dispersion are slightly 
lower for the various values of A in comparison to the 
unfrustrated case (Fig. and Fig. |3|3, left parts). 

The one-triplon spectral weight at a = 0.25 differs 
from the one at a = for momenta close to fc = tt 
(Fig. right part). The spectral weight is reduced for 
all momenta at a = 0.25 on increasing A. But the reduc- 
tion is smallest for k — In the limit of zero dimeriza- 
tion the one-triplon spectral weight Ii{k) vanishes for all 
momenta. The left part of Fig. |2t shows the one-triplon 
dispersion uj{k) for a = 0.5. The dispersion is highly 
reduced due to the increased locality of the triplon. At 
k = tt/2, the one-triplon state is an eigen-state of the 
system»2iS and it has an energy of J independent of A. 
In the limit of zero dimerization the system remains in a 
gapped stateiiSiSiiS^. 

The right part of Fig. O; shows Ii{k) for a ~ 0.5. 
The spectral weight is reduced for increasing A for all 
momenta except k = 7r/2. Since the one-triplon state at 
k = 7r/2 is an eigen-state independent of A its spectral 
weight is also constanii^. In contrast to the previous 
cases, there is also a one-triplon contribution for zero 
dimerization, at least for k = 7r/2, but most probably 
also in the vicinity of this momentum. 



B. Two- Triplon Contribution 

In this section we discuss the two-triplon contribution 
to the dynamical structure factor. The two-triplon con- 
tribution displays additional physics in comparison to the 
one-triplon part. The reason is that besides the kinetic 
part of the excitations also the triplon-triplon interac- 
tion is important and has to be included. An attractive 
interaction can lead to bound states of two triplons. Fur- 
thermore, the total momentum of two triplons does not 
fix the state of the system. There is also a relative mo- 
mentum between the triplons which is not fixed. Thus 
there is a continuum of two-triplon states for each given 
total momentum. Let us turn to the spectral properties 
of the two-triplon continuum and the two-triplon bound 
states. 

We decompose the action of the full observable 
^■^^^btrp in the two-triplon channel on the ground-state 
for fixed two-triplon momentum k in the two-triplon sec- 
tor 



i?^='|2trp|0) 



^Af=pi(fc,d)|fc,d). (12) 
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FIG. 3: One-triplon dispersion uj{k) (left panels) and one-triplon spectral weight Ii{k) (right panels) for A = {0.3; 0.4; 0.5; 0.6}. 
In Fig. (a) a = 0.0, in Fig. (b) a = 0.25 and in Fig. (c) a = 0.5. 



Here d denotes the relative distance between the two 
triplons and 

Al=;,ik, d) = -V2i J2 sin (^k{2l + l+d)y 

(13) 

The sum runs over all strong bonds I and is the 

amplitude for the creation of two triplons on dimers I 
and I + d hy Sf^. Here it is convenient to use a mixed 
representation in which the center-of-mass coordinate is 
Fourier transformed and the relative coordinate is dealt 
with in real space. The action of S]^ does not need to be 

R 5—1 L 5—1 

calculated due to the relation = —a_\_^ The 

basic unit length a/2 is again the distance between two 
neighboring sites. The spectra show results for momenta 
between and 1 in units of 7r/(a/2) (see Figl^J. In Figs. 
I4I7I the results for the two-triplon continua, the disper- 
sions of the two-triplon bound states and their spectral 
weights are shown. 

Fig. 2] displays the spectral density of the two-triplon 
continuum as a function of frequency and momentum. 
The spectrum for fixed momentum k is shifted by k in 
y-direction in order to provide a three dimensional view 
on the spectral densities. The lower and upper band edge 
is marked by solid gray lines. If there are any two-triplon 
bound states, their dispersion is also displayed as a black 
line. We denote the S = 1 two-triplon bound states as 
r„ and the S = two-triplon bound states as 5„ where 
n = {l,2,...}. 

Detailed information about the bound states is given 
in Fig. |S1 which consists of two parts for each parameter 
set. The left part shows an enlargement of the dispersion 
of the bound state and the lower bound of the two-triplon 
continuum. In the right part the corresponding spectral 
weight of the bound states is shown. The spectral weight 



is multiplied by the given factors for clarity. 

What are the general features of the obtained spectra? 
Due to the conservation of the total S'^-component there 
is no spectral weight at zero momentum. The energies of 
the system possess a reflection symmetry about fc = 7r/2 
which is a consequence of the inversion symmetry k <-!■ 
— fc and of the coupling of the momenta k and k + tt— . 
This symmetry can be clearly seen in the bound state 
energies and in the lower and the upper band edges of the 
continuum. It does not hold for the spectral weights'^^. 

In Fig.^Ji the spectral density for A — 0.3 and a = 0.0 
is shown. The spectral weight is mostly concentrated at 
the lower band edge of the continuum. There are two 
bound states centered about k = tt/2 leaving the con- 
tinuum at some finite momentum. The dispersions and 
the spectral weights of the bound states are plotted in 
Fig. E^. The points where the bound states are leaving 
the continuum can also be discerned in the singular be- 
havior of the lower band edge of the continuum. The 
spectral weight is mainly concentrated in the first bound 
state Ti . The spectral weight of the second bound state 
T2 is highly reduced. 

The binding energy of the bound states has its maxi- 
mum at fc = tt/2. It vanishes quadratically cx (fc — fcc)^ 
when the bound state enters the continuum. Correspond- 
ingly their spectral weight vanishes linearly cx | fc — fee | in 
accordance with the exemplary calculation provided in 
Ref.^^. There it was shown for square root type continua 
that the binding energy vanishes quadratically as func- 
tion of an external parameter while the spectral weight of 
the bound state vanishes linearly. The external parame- 
ter was the attraction strength; in the present case it is 
the total momentum which controls the relative strength 
of interaction and kinetic energy. 

Decreasing the dimerization to intermediate values 
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FIG. 4: Two-triplon spectral density l2{k,uj) for ^ with 
A = 0.3 and a = 0.0 (a), q = 0.25 (b) and a = 0.5 (c). Gray 
lines denote the lower and upper bound of the continuum. 
Black lines indicate dispersion of two-triplon bound states. 



(A = 0.6, see Figs. |HK and[7K), there are no qualita- 
tive changes in the spectrum. The spectral weight is 
mainly concentrated at low energies. The range of ener- 
gies, where the bound states exist, is slightly enhanced 
while the binding energy is slightly reduced. We expect 
that this tendency continues to lower values of the dimer- 
ization. For exactly zero dimerization, but not for an 
arbitrarily small but finite one, the first bound state Ti 
coincides with the lower bound of the continuum lead- 
ing to a square root divergence at the lower edge of the 
continuum for all momenta. Our expectation is strongly 
corroborated by the exact results for the spectral densi- 
ties in the sine- Gordon modelMiSi. The generic behavior 
is a square root behavior at the band edges. Only if 
a breather becomes degenerate with the multi-particle 
band edge the square root behavior switches to a square 
root divergence. Exactly the same characteristics ap- 
pears naturally in the triplonic description. Recall also 
that for the uniform spin chain at zero dimerization the 
square root divergence is well-known from the exact two- 
spinon contribution to the dynamical structure factos^. 

In Figs. 21d and ^ the spectra for strong dimeriza- 
tion and finite frustration a = 0.25 and a — 0.5 are 
shown. The corresponding information about the two- 
triplon bound states is plotted in Figs.Eb and|3;. The 
frustration makes the excitations more local and less dis- 
persive which leads to a narrowing of the two-triplon con- 
tinuum. At the same time the triplon-triplon interaction 
is enhanced causing an increase of the binding energy of 
the bound states. It can be nicely seen that for a — 0.25 
the first bound state extends over a wide range in mo- 
mentum space lying for small and large momenta very 
close to the lower band edge inside the continuum. For 
a = 0.5 the bound state Ti exists for all momenta. 

Due to the existence of the bound states only in a finite 
interval of momentum for both values a = and a — 
0.25 the qualitative distribution of their spectral weight 
is similar. This will be true for all values of frustration 
between and 0.25. The same holds for the bound state 
T2 for a = 0.5. In contrast, the spectral weight of the 
first bound state Ti, which is well separated from the 
continuum for a — 0.5, has its maximum at fc = tt. 

The whole two-triplon contribution, i.e. the two- 
triplon bound states and the two-triplon continuum, van- 
ishes for k — tt/2 a.t a — 0.5. Here the one-triplon excita- 
tion is an exact eigen-state of the spectrum and therefore 
comprises the total spectral weighlsSiS (see also preced- 
ing Section). 

We now turn to the influence of the frustration on the 
shape of the two-triplon continuum. In the case of vanish- 
ing frustration the spectral weight is distributed mainly 
close to the lower band edge for strong and intermediate 
dimerization. The spectral weight decreases monotoni- 
cally for higher energies. Turning on the frustration we 
observe a shift of spectral weight towards higher energies. 
In the case of strong dimerization this tendency is weak 
(Figs.^ and^) while for intermediate dimerization we 
observe a huge transfer of spectral weight (Figs.lHb and 



7 



1.00 



0.75 ■ 



0.50 - 



0.25 - 



0.00 



- 










^^;250 






^*_^50 




> 














(a) ■ 



1.8 1.9 
«[J] 



0.0 



\ ■ 


' ^^^^^ C.\J\J 














/ ■ 


(b) ■ 





1 ■ 1 ■ 




■ 


\ . 


/*25 












'•'"'12500 






^ 















0.5 
l(k) 



1.0 1.8 



1.9 
(0[J] 



0.0 



0.5 
l(k) 



1.0 



1.85 1.90 
(0[J] 



0.0 



0.5 
l(k) 



1.0 



FIG. 5: Two-triplon bound states for R^^^ with A = 0.3 and a — 0.0 (a), a — 0.25 (b) and a — 0.5 (c). Left panels show the 
dispersion of the bound states; right panels show the spectral weights of the bound states multiplied by the indicated factors. 
Gray lines denote lower bound of the continuum. 



Et). This transfer produces a non-monotonic shape for 
intermediate dimerization and a = 0.25 having a mini- 
mum of spectral weight inside the continuum (Fig. |B|d). 
Increasing the frustration {a = 0.5) shifts the minimum 
to the lower band edge. The spectral weight is mainly at 
the upper band edge (Fig.int). 

These observations are very similar to the results ob- 
tained by exact diagonalization at finite temperatures for 
the dynamical structure factor of a homogeneous spin 
chain including frustration^i. There a decrease of spec- 
tral weight at A; = TT inside the continuum is observed on 
increasing the frustration. This results in a high-energy 
maximum for large frustration and a minimum inside the 
continuum, i.e., a trough- like shape. From this compar- 
ison we conclude that our findings represent the generic 
features which are also valid in the limit of vanishing 
dimerization. 



C. Comparison with Field Theory Results 

In the first part of this section we concentrate on the 
S — 1 breather and the corresponding behavior of the 
lower band edge of the two-triplon continuum. In the 
second part we have a closer look at the energy of the 
S — breather and at the one-triplon gap. Finally, the 
importance of marginal terms for the quantitative shape 
of spectral densities is discussed. A detailed analysis of 
the S ~ channel which contains also a discussion of the 
singlet two-triplon bound states is presented in the next 
section. 

Let us first look at the results obtained from bosoniza- 
tion and the continuum limit renormalization group 
approachSS*^. The bosonized form of the dimerized and 



frustrated spin chain reads 



FT 



V 

2^ 



K {TTllf + K-^ {d^<pf dx 

/OO 
[Mcos(2$) + Dcos(4$)]da;, (14) 
-OO 



where v is the spin- wave velocity and K — 1/2 the in- 
teraction parameter for the isotropic chain. The value 
of the critical frustration q;o,c = 0.241167 depends on 
the physics at short distances and is only accessible by 
numerical techniquesi^iSiSiiS^. The i5 cos(2<I>) term is 
strongly relevant while the D cos (4$) term is marginally 
irrelevant for ao < ao,c and marginally relevant for 
ao > q;o,c- It is commonly accepted and numerically 
confirmed^SiSLSiS that the marginal term can be ne- 
glected best for ao = ao.c (for further discussion see be- 
low). Then a simple sine-Gordon model at K = 1/2 
remains of which the spectral densities are knownS^iS. 

The 5" = 1 response function displays a square root di- 
vergence oc {io — ujQ)~^/^ at the lower band edge. Here the 
energy of the S — I breather ix'br,s=i is exactly degen- 
erate with the lower band edge {uJhi-,s=i = 2A). This is 
in agreement with what we find at a = 0.25, cf. Figs.^jD 
and I^Id. Without any frustration, however, we find a 
square root behavior cx (w — tooY^'^ , cf. Figs.0^ andE^. 
Hence we conclude that the sine- Gordon model does not 
describe the unfrustrated, dimerized spin chain exhaus- 
tively. But the sine-Gordon model applies to the spin 
chain at critical frustration where the higher cosine-term 
cos(4$) vanishes. It is interesting to note that we find a 
square root divergence for a « ao,c and not for ao ~ ao.c- 

The conclusion about the applicability of the sine- 
Gordon model to the unfrustrated and dimerized spin 
chain is in agreement with the results of the numerical 
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FIG. 6: Two-triplon spectral density l2{k,uj) for ^ with 
A = 0.6 and a = 0.0 (a), q = 0.25 (b) and a = 0.5 (c). Gray 
lines denote lower and upper bound of the continuum. Black 
lines indicate dispersion of two-triplon bound states. 



investigation of the bound states^. It is known that 
the elementary excitations of the SU(2) symmetric sine- 
Gordon model consist of soliton and antisoliton excita- 
tions and two breathers, bound states, plus one breather 
which is degenerate with the lower band edgei*^2iS£. The 
lowest-lying breather is degenerate with the soUton and 
antisoUton excitations and corresponds to the 5*^=0 
triplet state in spin language. This fixes the interaction 
parameter to = 1/2. The second breather is assigned 
to a spin singlet excitation, since there is no counter part 
in the soliton or antisoliton sector. The ratio between the 
energy of the 5 = breather a^br,s=o and the one-triplon 
gap A is exactly a/3 a,t K — 1/2. 

In Fig. |H1 this ratio is shown for various values of 
A — {0.3; 0.4; 0.5; 0.6} versus the bare frustration aQ. We 
have used Fade extrapolants for A = 0.6. It can be clearly 
seen that almost all points (except the case for strong 
dimerization A = 0.3) fall onto one curve. The values 
for ofg c and are included as solid lines. The point 
where these two lines cross lies on the calculated curve 
and corresponds to the prediction of the SU(2) symmetric 
sine-Gordon model. In all other cases (oq ^ q;o,c) the ra- 
tio ajbr,s=o/A differs from \/3. This is due to corrections 
resulting from the marginal term cos(4(E>). Our calcula- 
tions agree perfectly with previous numerical results'^^. 

The relative importance of the two cosine-terms in 
Eq. H14|) are at finite dimerization is a subtle issue which 
we discuss in the following. 



In the self-consistent harmonic approximatio] 



,61.62 



we 



replace $ <&ciass + ^fluct where only the fluctuation 
part is of operator character. The cosine-terms can then 
be approximated by 



(15) 



cos(n<I>) 

exp{-{n^/2)a{x))cos{n^cUss{x)) [ 1 - y^fl^ct 



where a{x) := ($fluct(2^)^)- This kind of approach corre- 
sponds to renormalization in first order. In the ground 
state (without solitons) one has $ciass — 0. In the un- 
gappcd phase the fluctuations diverge, hence a ~* oo. 
But in the gapped phase, the fluctuations are cutoff at 
low energies so that a = — (iir/2) ln(A/Ao) where A is 
the gap and Aq is proportional to the ultraviolet cutoff. 

Since the square of the gap A^ is proportional to the 
coefficient of $(x)^ one obtains from the term (5cos(2$) 
the self-consistency equation 



A^ cx (5exp(-2cr) 
A cx <5i/(2-K) 



(16) 



yielding the well-known A cx (5^/"^ by Cross and FisheisS^. 
Hence the total contribution of this cosine-term is pro- 
portional to A^ or 5'^^^ . The crucial point to note is that 
the amplitude of the second cosine-term cos(4$) is of the 
same magnitude exp(— Scr) cx A'^^ which also yields A^ 
or (5^/^ for K = 1/3S^. Hence the self-consistent har- 
monic approximation tells us that even in the regime 
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FIG. 7: Two-triplon bound states for R^^^ with A = 0.6 and a — 0.0 (a), a = 0.25 (b) and a = 0.5 (c). Left panels show the 
dispersion of the bound states; right panels the spectral weights of the bound states multiplied by the indicated factors. Gray 
lines denote lower bound of the continuum. 
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where the frustration is marginally irrelevant a < ac 
it influences the low-energy physics on the quantitative 
level. 

Considering, however, the renormalization to second 
order the marginally term is reduced logarithmically, 
hence the name "marginally irrelevant" (Refs. 34,64). 
The flow is valid down to the infrared cutoff which is in 
our case the energy gap due to the dimerization. Note 
that we are working in the thermodynamic limit L = oo. 



Hence the suppression of the marginally irrelevant cosine 
term to the relevant one is logarithmically in 6 and so 
very slow. Thus it is possible that the scales which we 
discuss in this paper {6 down to 0.25) are still too large 
to see the emergence of the physics of a pure sine- Gordon 
model with one cosine only. But we find it puzzling that 
no precursors of the convergence to the predictions of 
the pure sine-Gordon model can be seen. The ratio of 
the energy gaps shown in Fig. |H| (or those obtained in 
Ref. 37) appears almost idependcnt of S as far as we 
could investigate. 

Note that the vanishing of the cos(4$) term at the end 
of the flow of / is not identical to the question whether this 
term is present at the beginning of the flow. We presume 
that this difference explains the puzzling finding that the 
bare coefficient of the Umklapp term cos(4$) vanishes^ 
at ao — 1/6, that means relatively far from the quantum 
critical point ao = 0.241. 

Our results in Figs. I4l6r a-b') show that the square 
root divergence known from the sine-Gordon model^^*^ 
is changed to normal square root behavior. So the quan- 
titative changes of the low-energy Hamiltonian influence 
the shape of the spectral densities qualitatively. Hence 
for spectral densities one must know whether the effective 
low-energy is a (simple) sine-Gordon model or a double 
sine-Gordon model with two cosine terms. The answer 
depends on the energy scale considered, i.e. the value of 
S. 

For strong frustration a = 0.5, cf. Figs.^J; and|Sl:, the 
physics is dominated by bound states. Their number pro- 
liferates for decreasing dimerizatior&S&. At 5 = there 
are inflnitely many bound states densely distributed be- 
tween A and 2A. They form the continuum which can 
be understood as two-spinon continuum22iS4S&. The val- 
ues of A for which we display the spectral densities in 
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Figs. 21: andEJ; are still too low, i.e., too far in the dimer- 
ized regime to see more than two bound states. This was 
also observed by Zheng et al.^"^ . They found in a series 
expansion up to order A^^ three singlet and three triplet 
bound states. 

We attribute the fact, that only a small number of 
bound states could be found so far, to the limited range of 
the effective interaction. In the perturbative approaches 
for a = 0.5 an order of A" corresponds to a maximum 
range of [n/2]. So even calculations at n = 19 provide 
only a potential of small finite range which does not allow 
for many bound states. The alternative presumptionSi 
that the lacking bound states is found in the channels 
with more than two triplons would require that the spec- 
tral weight is passed to channels with more and more 
triplons. No channel with only a finite number of bound 
states may retain a finite spectral weight at zero dimer- 
ization since at zero dimerization only a continuum is 
found^^. None of our results is in favor of this scenario 
so that we are convinced that the range of the interaction 
is the crucial point. The fact that an expansion to higher 
order finds more bound states in the two-triplon sector 
supports the view that the range of the interaction mat- 
ters. But the precise description of the deconfinement 
transition for vanishing dimerization is still an open is- 
sue. Future developments like self-similar realizations of 
the continuous unitary transformations will help to re- 
lease this constraint on the range of the interaction^SiS&. 

Concerning the present data at strong frustration, an 
important trend is that on increasing A the spectral 
weight is shifted towards higher energies. 



V. S=0 EXCITATIONS 

In this section we concentrate on the two-triplon con- 
tribution to the spectral density with total spin zero 
which is relevant for optical experiments. For 5*^0 
one needs at least two triplons which form together a 
state with vanishing total spin. For the case of isolated 
dimers (A = 0) the total spectral weight is in the two- 
triplon contribution. At finite A the spectral weight in 
the two-triplon channel is reduced and spectral weight is 
also found in channels with more than two triplons. 

In the limit of vanishing dimerization for zero frus- 
tration we can show in a similar analysis as for the dy- 
namical structure factor that the two-triplon contribu- 
tion possesses almost the total spectral weigh1i224&. The 
two-triplon contribution is therefore the only sizable con- 
tribution for the whole range of dimerizations for the un- 
frustrated case. In presence of frustration the analogous 
analysis is quantitatively more difficult as stated before. 
But there are again indications that only a small number 
of triplons dominates the spectral properties. Therefore, 
we investigate the leading two-triplon contribution for 
the 5 = channel. For the dimerizations treated in this 
work there is no doubt that the two-triplon contribution 
is the only sizable term. The crucial point, however, is 



to which extent we can obtain the generic features which 
govern also the limit of vanishing dimerization. All re- 
sults obtained so far show that the one- and two-triplon 
contributions capture indeed the relevant physics. 

The physical observable R^^^ for S = excitations 
locally reads 

i.e., it is a sum over nearest neighbor (NN) and next 
nearest neighbor coupling (NNN). The coefficient /3 is a 
measure for the relative different strength between the 
two couplings. It depends on the underlying microscopic 
physics and will not be discussed in this work. As illus- 
trated in Fig. 1^ these observables are given by 

Rio=%N = (1+7)So,lSo.r + (1-7)So,rSi,l (18) 
for nearest neighbor (NN) coupling and 

Ao7,NNN — SoxSi,L + So,rSi,r (19) 

for next-nearest neighbor (NNN) coupling where 7 is pro- 
portional to the dimerization 6. -Rjo^nn ^ '^'^^^ 
couplings on weak and strong bonds. 

We will restrict our discussion to the case of near- 
est neighbor coupling on the weak bonds and on the 
next-nearest neighbor coupling. This choice is motivated 
by the relevance of various observables for Raman spec- 
troscopy and infrared absorption in the limit of vanishing 
dimerization. Raman spectroscopy measures excitations 
with total momentum zero while infrared absorption is 
governed by the response at large momenta^^iZS. 

In the case of a uniform spin chain without frustration 
the nearest neighbor Raman operator commutes with 
the Hamiltonian and one obtains a vanishing Raman re- 
sponse. Therefore the next-nearest neighbor Raman op- 
erator is the leading contribution in terms of a Loudon- 
Fleury scattering theor y'^^i'''^ . In contrast i?NN does not 
commute for finite momenta with the Hamiltonian and 
will be the most important contribution to the infrared 
absorption. For simplicity, we do not treat -Rnn com- 
pletely but only the weak-bond part. This is no major 
restriction because we are interested in the generic prop- 
erties of these quantities. In addition, the weak-bond 
part dominates for strong dimerization. 

We will discuss these two observables i?NN and -Rnnn 
separately. For a direct comparison with experimental 
data one should take the sum over all contributing parts 
of i?'^^°(fc) to account for possible interference effects. 
The necessary superposition, however, depends strongly 
on the details of the system and cannot be discussed gen- 
erally. 

The action of the full observable on the ground-state 
is decomposed again for fixed total momentum k in the 
two-triplon sector by 

^NN°woak(^)|0) = ^ ^2trp,NN,wcak(''^; '^)l^i (20a) 
d 

^NNN(fc)lO) = Y.^2tt,N^^ik,d)\k,d) , (20b) 

d 
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where 



A. Case R'^^° 



^ft7p,NN,weak(fc,rf) = V2 J] a-^f/'' COS (^(2/ + d) ) (21a) 

^2trp,NNN(^i = (21b) 

J2V2allf^^cos{k{2l + l/2 + d)) . 



Here d is the distance between the two triplons, a^i+d '^^ 



-1/2 I 



1/2 1 3/2 

' ^NNN 



— a — ^ 



FIG. 9: Sketch of the local observables for S—O excitations. 
-RnJ'^ is a sum of couplings on strong bonds (double lines) and 
weak bonds (single line). -Rnnn couples next nearest neighbor 
spins. 



is the amplitude for the creation of two triplons on the 
dimers I and I + dhy So,rSi_l, and aj'j!^^^ is the ampli- 
tude for the creation of two triplons on the dimers I and 
l + dhy So,lSi,l- It is not necessary to calculate aj^;'^^^ 



due to the symmetry aj';^^^ 



a 



L,NNN 
-l- 



-V 



The basic unit 



length a/2 is again the distance between two neighboring 
sites. Momentum is measured in units of 7r/(a/2). 

We first discuss the symmetries of the two observables. 
-^NN^wcak possesses a reflection symmetry about k = n/2 
for the same reasons as the Hamiltonian. Any mode at 
k which is created by ^NN°woak coupled to the mode 
at k + TT. In addition every mode at k corresponds to 
the mode at —k. Therefore the whole spectral density 
will be symmetric about k = 7r/2. This symmetry is 
absent for i?NNN- ^or i?NNN *be spectral weight is mainly 
concentrated at small and intermediate momenta while 
it vanishes exactly at fc = tt. The latter follows from 
the fact that at fc = tt the observable creates an odd 
state with respect to reflection about the axis 1/2 (see 
Fig. Owhile a singlet made from two triplets is always 
an even state with respect to particle exchange. 

In Figs. llOlT^ the spectral densities for -RNN^wcak ^^'^ 
shown and in Figs. [TTll7l the corresponding densities for 
Let us first discuss the results for the nearest 
neighbor (NN) coupling ^NN'weak passing then to the re- 
sults for next-nearest neighbor coupling -R^nn- Finally, 
the implications for Raman spectroscopy and infrared ab- 
sorption will be assessed. 



In Fig. llUb the spectral density for strong dimerization 
(A = 0.3) and vanishing frustration is depicted. The 
corresponding information about the dispersion and the 
spectral weight of the two-triplon bound states is shown 
in Fig. Illb . The same notation as in the Section for 
S = 1 excitations is used. The spectrum is symmetric 
about k — 7r/2 due to the inversion symmetry k ^ —k 



'^) = figN°wcak(^)- 



and i?|N°wcak(fc 

We find two S — two-triplon bound states 5*1 and 5*2. 
The triplon-triplon interaction is larger in the total 5 = 
channel than it was in the total S = 1 channel. There- 
fore the binding energy of the bound states is enhanced 
and the first bound state Si exists for all momenta in 
contrast to the 5=1 case. In general, the S = channel 
is dominated by the bound states which carry most of 
the spectral weight. This statement applies also to the 
experimental relevance, see below. 

The dispersion of the two-triplon bound state Si is 
roughly sinusoidal having three extrema at momenta fc = 
{0;7r/2;7r}. The binding energy is largest for k — Tr/2 
while it becomes small near momentum zero and tt. The 
spectral weight of is roughly proportional to the bind- 
ing energy. The second singlet two-triplon bound state 
S2 exists only in a finite interval about k = tt/2. The 
spectral weight of S2 vanishes at k = tt/2 and possesses 
two maxima below and above k = 7r/2. 

The spectral weight of the two-triplon continuum is 
concentrated at small frequencies. At small and large 
momenta this effect is enhanced due to the vicinity of 5*1 . 
Lowering the dimerization we see no qualitative changes 
in comparison to the case of strong dimerization fFig.ll2b. 
and Fig. Il^b). So we expect that the dispersion of the 
bound state Si is degenerate with the lower band edge of 
the two-triplon continuum inducing a square root diver- 
gence. This expectation is supported also by numerical 
results for the second breather^^I. 

For clarity, we like to emphasize again that one must 
clearly distinguish the the case of zero dimerization and 
the case of small, but finite, dimerization. For zero dimer- 
ization bosonization predicts a, 1/uj divergence at fc = tt 
which becomes {lu — cjo)^^^^ close to fc = TiiS. This 
has been used for instance in the empirical calculation of 
Lorenzana and Eder— . For small, but finite, dimeriza- 
tion the sine-Gordon model prediction of a square root 
behavior without divergence applies to the critical frus- 
tration and in the region around this value. For other 
values of the frustration a breather may coincide with 
the lower band edge implying a square root divergence. 

In Figs. llOb-c the spectral density at finite frustration 
for strong dimerization is shown. At a = 0.25 we find 
three bound states Si , S2 and S3 . The binding energy of 
Si increases drastically on turning on the frustration, es- 
pecially at small and large momenta. The spectral weight 
increases in a similar fashion for these momenta. The 
third two-triplon bound state 6*3 exists merely in a very 
small region about fc = 7r/2. The binding energy and 
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the spectral weight are tiny. The spectral weight has a 
maximum at fc 7r/2. 

The two-triplon continuum for a — 0.25 does not show 
much structure. This is a consequence of the fact that at 
almost all momenta no bound state is close to the lower 
band edge of the continuum. Decreasing the dimerization 
no qualitative changes are seen (see Fig. 112b). 

At a = 0.5 we detect two bound states. The dispersion 
of Si becomes flatter which again holds also for the spec- 
tral weight distribution. The biggest change can be seen 
in ^2 . This bound state exists for almost all momenta in 
contrast to the cases a — {0;0.25}. In the regions close 
to k — and symmetrically close to fc = tt, the bound 
state ^2 does not exist, but it can be thought to lie just 
above the lower band edge implying an almost divergent 
behavior of the two-triplon continuum. 

Smirno\i2^i^ showed that the corresponding spectral 
density of the sine-Gordon model displays a square root 
behavior at the lower band edge. This applies to the 
S — channel of the frustrated spin chain at a ^ ao,c 
if the marginal term D cos(4$) of Eq. (|14() is neglected. 
This neglect is quantitatively justified for a = ag.c- In- 
deed, our results clearly show a square root behavior for 
a = 0.25. As for the 5 = 1 case, we find that the pre- 
dictions of the sine-Gordon model for the physics of the 
spin chain are verified for a — ao.c- For other values, 
notably a ~ and a = 0.5, we find square root diver- 
gences or strong tendencies towards square root diver- 
gences. Again, such divergent behavior results from the 
vicinity of bound states, here in the S = sector. 



B. Case 7?S^n 




CM 

is 
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In Figs.lT^andlTHIthe spectral densities of the observ- 
able i?NNN various values of dimerization and frus- 
tration are shown. The information on the singlet two- 
triplon bound states is plotted in Figs. El and El AH 
considerations concerning the energetic properties of the 
S = channel are the same as for -R^nn ^^'^ need not 
be discussed again in this section. 

We concentrate on the spectral differences between 
the two observables. The spectral weight is suppressed 
for large momenta due to symmetry reasons. This can 
be clearly seen for the two-triplon continuum and the 
spectral weight of the two-triplon bound states. There- 
fore only momenta k £ [0, 7r/2] carry significant spectral 
weight. In this region we do not find qualitative changes 
to the results for of ^NN^wcak- 

The most important difference is a change of the spec- 
tral weight distribution of 5*1. At a = the spectral 
weight has a maximum at fc = 7r/2 which is similar to the 
case of i?NN'wcak- Finite frustration shifts the maximum 
to fc = 0. Close to the critical frustration the spectral 
weight is almost constant for k S [0, 7r/2]. At a = 0.5 
the spectral weight is a monotonically decreasing func- 
tion (from fc = to fc = tt). 




Q)[J] 



FIG. 10: Two-triplon spectral density l2{k,uj) for -RNN^wcak 
with A = 0.3 and a = 0.0 (a), q = 0.25 (b) and a = 0.5 (c). 
Gray lines denote lower and upper bound of the continuum. 
Black lines indicate dispersion of two-triplon bound states. 
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FIG. 11: Two-triplon bound states for -RNN^woak with A = 0.3 and a — 0.0 (a), a = 0.25 (b) and a = 0.5 (c). Left panels show 
the dispersion of the bound states; right panels the spectral weights of the bound states multiplied by the indicated factors. 
Gray lines denote lower bound of the continuum. 



C. Raman Spectroscopy 



D. IR- Absorption 



The dominant observable for magnetic light scatter- 
ing (Raman response) using the standard Fleury-Loudon 
scattering theorjj£i*£^ is 



E 



R 



s=o 

loc,NN 



{i) + I3R 



■locNNN 



(0) 



(22) 



where the sum runs over all spins. The Raman response 
is therefore the A: = contribution to the spectral density 
we have discussed in the last section. We focus here on 
the case of next-nearest neighbor coupling which is the 
leading process in the case of a uniform Heisenberg chain 
without frustration. 

In Fig. ^1 the Raman response for next-nearest cou- 
pling is shown at zero frustration (a), close to critical 
frustration a = 0.25 (b) and for a = 0.5 (c). In each 
graph the spectrum is shown for A — {0.3; 0.4; 0.5; 0.6}. 
In these figures, a broadening of F = 0.01 is used and the 
spectra are shifted in y-direction for clarity. The spectral 
densities for a = are multiplied by 6. 

For the dimerizations considered here, the spectra are 
dominated by the first S — Q two-triplon bound state 
Si. This dominance is enhanced by the frustration. In 
Fig. I18h , the case of vanishing frustration is shown. Due 
to the finite broadening and the small binding energy of 
the bound state there is no separation of the two-triplon 
bound state and the two-triplon continuum. An increase 
of A reduces the weight of the bound state and gives rise 
to a broad featureless continuum. 

In Fig. 118b and 118b the results for finite frustration 
are plotted. The binding energy of the bound state 5*1 is 
enhanced and one can clearly separate the contribution 
of the bound state 5*1 and the continuum. The spectral 
weight of the two-triplon continuum is very small. 



In this section we apply our results to phonon-assisted 
infrared absorption of magnetic excitationaSSiZS. This 
technique allows to study the spin-spin correlation func- 
tion by measuring the optical conductivity. The direct 
absorption of two magnetic excitations is generically not 
allowed due to inversion symmetry. However, this se- 
lection rule can be broken by simultaneously exciting a 
phonon. The leading infrared-active magnetic absorp- 
tion is a two-triplon-plus-phonon processS^iZS. Due to the 
momentum of the excited phonon, the magnetic spectra 
I{k,Lu) have to be integrated over all momenta weighted 
with a phonon-specific form factor. 

The absorption spectra are sensitive to the S = two- 
triplon bound states. Especially the extrema yield promi- 
nent van-Hove singularities in the density of states which 
can be identified in experiment. In this way, the first 
experimental evidence for the two-triplon bound state in 
cuprate spin ladder systems2^ was possible. We therefore 
expect interesting line shapes in the optical absorption 
also for dimerized and frustrated spin chain systems. 

The absorption coefficient is given b}i24 



a{Lu) — aoUjI^^Xuj — ujq) 



(23) 



Here oq is a constant depending on the material and ujq 
is the phonon frequency. The phonon is considered to be 
local and without dispersion. The function /^^ is given 
by 



r™(w) = 167r^sin4(fc/2)/(/c,w). 



(24) 



The specific form factor given is strictly valid only for 
a uniform Heisenberg chain. It was successfully used to 
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FIG. 12: Two-triplon spectral density l2{k,uj) for -RNN^wcak 
with A = 0.6 and a = 0.0 (a), a = 0.25 (b) and a = 0.5 (c). 
Gray lines denote lower and upper bound of the continuum. 
Black lines indicate dispersion of two-triplon bound states. 



explain the optical absorption in uniform cuprate spin 
chains^^iS^. We use the same form factor also for the 
dimerized and frustrated chain in order to explore the 
general features of the optical conductivity and to com- 
pare the line shapes at finite dimerization with the line 
shapes at zero dimerization. In a detailed analysis of 
experimental data one must analyze which phonons are 
involved and which specific form factors matter. 

In Fig. I19h,-c and Fig. I2()b-c the optical absorption 
a{uj) for various dimerization and frustration is shown 
for i?NNN ^i^d ^NN^wcak- Here oo is set to one and wo 
is set to zero. The spectra are plotted with a broad- 
ening of F = 0.01 which is a reasonable value in view 
of experimental resolutions. In the insets we show the 
contribution of the two-triplon continuum without the 
broadening to highlight the shape of the continuum con- 
tributions and to distinguish it from contributions of the 
bound states. 

The phonon form factor favors large momenta while it 
reduces the contribution of small momenta. Hence the 
discussion of the spectral densities implies that ^NN*wcak 
is more relevant than i?NNN ^^^r a{uj). This can also be 
seen in the absolute heights of the spectra in Fig. 1 191 and 
Fig. [501 In addition, we expect that the nearest neighbor 
coupling is stronger than the next-nearest neighbor one 
because exchange processes of longer range will generi- 
cally be less important. 

We start our discussion with i?NN'wcak- stated 
above the first S — two-triplon bound state Si car- 
ries most of the spectral weight for all momenta. Hence, 
it is of crucial importance for the optical absorption. In 
Figs. [TTl and [T^ we show that the dispersion Wbound(^) 
of the bound state 5*1 possesses three extrema at A: = 
{0, 7r/2, tt}. So we obtain three van-Hove singularities in 
I^^ {uj). The spectral density is symmetric about k = it/2 
so that two van-Hove singularities coincide and there are 
two peaks resulting from the bound state Si in the opti- 
cal conductivity. The weight of the minimum at A: = is 
suppressed by the phonon form factor. This implies that 
the regions about k — t: /2 and k = tt dominate. 

In Fig. 119b the optical absorption for a dimerized chain 
(A = {0.3,0.3,0.5,0.6}) without frustration is depicted. 
The spectra are shifted in y-direction for clarity. The 
line shape is dominated by a small peak at low energies, 
a sharp peak at intermediate energies and a broad struc- 
ture at high energies. The first two features are mainly 
produced by the above mentioned van-Hove singularities 
resulting from the extrema of the bound state dispersion 
of 5*1 . The second peak is dominant because the spectral 
weight has a maximum for k = tt/2, see Fig.lllb andll3h,. 
For increasing A this peak looses intensity while the first 
peak becomes more pronounced. The latter effect is due 
to the increasing binding energy of S*! at k — tt. 

For strong dimerization the feature at low energies, 
which is more like a shoulder than like peak, is an ef- 
fect of the two-triplon continuum (inset Fig. I19bl. The 
second bound state S2 is of no greater relevance for the 
optical absorption because it has zero spectral weight for 
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FIG. 13: Two-triplon bound states for -RNN^wcak with A = 0.6 and a — 0.0 (a), a. — 0.25 (b) and q — 0.5 (c). Left panels show 
the dispersion of the bound states; right panels the spectral weights of the bound states multiplied by the indicated factors. Gray 
lines denote lower bound of the continuum. 



k — IT /2 which is the only extremum of the bound state 
dispersion. In addition, the binding energy is very small 
without frustration and so is the corresponding spectral 
weight. 

Lorenzana and Eder— calculated the two-spinon-plus- 
phonon contribution to the optical absorption for a uni- 
form Heisenberg chain. The line shape consists mainly 
of three parts: a concave uprise at low energies which 
vanishes for zero frequency, a singularity at intermediate 
energies and a convex tail for higher frequencies. It is 
very interesting to see that all these features have pre- 
cursors at finite dimerization which are captured in the 
triplon picture. 

In the limit of vanishing dimerization the system be- 
comes gapless and the spectra therefore start at zero en- 
ergy. As long as there is some finite dimerization the 
bound state Si exists and produces the concave uprise at 
small energies and the singularity at intermediate ener- 
gies resulting from the maximum at fc = 7r/2. We expect 
that for vanishing dimerization (A = 1) the dispersion of 
Si coincides with the lower band edge of the two-triplon 
continuum leading to a square root divergence at the 
lower band edge for all momenta. Since the dispersions 
and the band edges display an extremum at k — it /2 this 
divergence leads to the singularity discernible at interme- 
diate energies. The convex tail at the upper band edge 
is equally present even for strongly dimerized chains, see 
inset in Fig. 119b. It is a consequence of the convex square 
root behavior at the upper edge of the two-triplon con- 
tinuum. 

In Figs. 119b and 119b the optical absorption at fi- 
nite frustration a — 0.25 and a — 0.5 for the same 
value of dimerization is shown. As discussed earlier the 
frustration enhances the triplon-triplon interaction and 
increases the binding energy of the two-triplon bound 



states. As can be clearly seen in Figs. Illb-c and 
Figs. I13b-c the spectral weight of 5*1 at /c = tt increases 
compared to the weight at A: = 7r/2. Therefore, the first 
peak in the optical absorption becomes more and more 
prominent on increasing frustration. This leads to the 
most important features at large frustration. At a = 0.5 
the spectral weight of 5*2 is also sizable. Besides the con- 
tribution of the two-triplon continuum (inset of Fig. 119b') 
an additional peak appearing for decreasing dimerization 
can be discerned. 

The optical absorption for i?NNN plotted in Fig. 1201 
In the insets an enlargement of the line shapes is depicted 
in order to highlight fine structures. The main difference 
to the discussion of the optical absorption produced by 
-RNN,wcak are the consequences of the different symme- 
tries of the observables. i?NNN suppresses the spectral 
weight for large momenta. Thus the optical response is 
weak due to the phonon form factor which stresses large 
momenta. In addition, the van-Hove singularity result- 
ing from fc = TT of 5*1 is suppressed so that only a weak 
shoulder can be observed at low energies, independent of 
frustration and dimerization. At a = 0.25 the additional 
side structures are produced by the bound states S2 and 
^3. 



VI. SUMMARY 

In this work we have presented results for the spec- 
tral densities of the dimerized and frustrated Heisenberg 
chain. We used a perturbative realization of the con- 
tinuous unitary transformations starting from the limit 
of isolated dimers. By the transformations an effective 
model is obtained which conserves the number of triplons. 
Results for the one-triplon and the two-triplon contribu- 
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FIG. 14: Two-triplon spectral density l2{k,uj) for ii^NN with 
A = 0.3 and a = 0.0 (a), a = 0.25 (b) and a = 0.5 (c). Gray 
lines denote lower and upper bound of the continuum. Black 
lines indicate dispersion of two-triplon bound states. 



tion to the spectral density have been shown for strong 
(A — 0.3) and intermediate (A = 0.6) dimerization and 
for various values of the frustration {a = 0; 0.25; 0.5) and 
for total spin one and zero. 

In the first part of this paper we examined the dy- 
namical structure factor which is relevant for inelastic 
neutron scattering experiments. The one-triplon contri- 
bution contains most of the spectral weight at strong 
and intermediate dimerization. We provided results for 
the one-triplon dispersion uj{k) and the fc-resolved spec- 
tral weight Ii{k). The one-triplon dispersion becomes 
larger on lowering the dimerization while it becomes flat- 
ter on increasing the frustration. The spectral weight 
Ii{k) is mainly concentrated at fc = tt. In the limit A — > 1 
the one-triplon contribution vanishes except for a = 0.5 
around k = tt/2. 

Subsequently we discussed the two-triplon contribu- 
tion to the dynamical structure factor. We have provided 
results for the spectral density of the two-triplon contin- 
uum and for the dispersion and the spectral weight of the 
two-triplon bound states. 

For the unfrustrated spin chain, the spectral weight 
is concentrated at the lower band edge at larger mo- 
menta. Two-triplon bound states only exist in a finite 
region about k — 7r/2. Increasing the frustration leads to 
a shift of spectral weight to higher energies. At A = 0.6 
and a = 0.5, the spectral weight is shifted almost totally 
to the upper band edge at k — n. This transfer of spectral 
weight is also found for A = 1 by exact diagonalization 
at finite temperatures^ . 

The behavior of the lower band edge changes strongly 
by the variation of the frustration. Generically, we find 
a square root behavior of the lower band edge. It is 
a consequence of the hardcore interaction between the 
triplons which makes it difficult for them to pass each 
other. In contrast to this finding, wc find a square root 
divergence at the lower band edge for a = 0.25. Here the 
energy of a two-triplon bound state is degenerate (to the 
precision of our analysis) with the lower band edge of the 
two-triplon continuum. 

We compared the latter finding with results obtained 
from field theory'^^. In contrast to our finding for the 
unfrustrated spin chain, field theory predicts a square 
root divergence for the lower band edge of the dynamical 
structure factor of the sine-Gordon model. Thus the com- 
monly used reduction of the spin chain to a sine-Gordon 
by neglecting the marginal operator cannot be justified 
quantitatively, at least not for the values of dimerization 
considered here. We showed that in the self-consistent 
harmonic approximation (corresponding to renormaliza- 
tion in first order) the marginal operator D cos (4$) is as 
important as the mass operator (5cos(2$) at any finite 
dimerization. Renormalization in second order, however, 
predicts a slow logarithmic suppression of the marginal 
term24 below the critical frustration. 

We find that square root behavior represents the 
generic behavior. A square root divergence occurs if 
a two-triplon bound state is degenerate with the lower 
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FIG. 15: Two-triplon bound states for -Rnnn with A = 0.3 and a — 0.0 (a), a = 0.25 (b) and a = 0.5 (c). Left panels show the 
dispersion of the bound states; right panels the spectral weights of the bound states multiplied by the indicated factors. Gray 
lines denote lower bound of the continuum. 



band edge of the continuum. In the field theoretic lan- 
guage this degenerate bound state, which has not yet 
emerged from the continuum, is the third breather. We 
find the concomitant square root divergence for a w a^ c- 

The applicability of the sinc-Gordon model to the un- 
frustrated, but dimerized spin chain is further questioned 
by the study of the excitation energies of the bound 
states. The ratio of the excitation energies of the S = 
two-triplon bound state and the one-triplon gap is ex- 
actly ^/i in the SU(2) -symmetric sine-Gordon model. 
We find this ratio only for ao = ao,c in agreement with a 
previous numerical studjiSi. At present, we do not know 
why ao = Q^o,c is required to retrieve the field theory re- 
sult for the second breather, but a « ao,c to retrieve the 
field theory result for the third breather. We cannot ex- 
clude completely that one has to go to much lower values 
of the dimerization, i.e., closer to A = 1. But in view of 
the fact that the marginal operator cannot be neglected 
for a quantitative description, we presume that the ef- 
fective low-energy model to be considered is the double 
sine-Gordon model for which the SU(2) symmetry con- 
dition will be different from K ~\/2 and hence also the 
ratios of the breather energies will differ from the ratios 
known for the sine-Gordon model^*^. 

In the second part of this work, we discussed spectral 
properties of the dimerized and frustrated spin chain for 
excitations with total spin zero which are relevant for 
optical experiments. We presented results for the two- 
triplon contribution which contains most of the spectral 
weight. 

We examined two different observables: a nearest 
neighbor (NN) coupling on the weak bonds i?NN°woak ^"^^ 
a next-nearest neighbor (NNN) coupling -R^nn- The ob- 
servables obey different symmetries. The nearest neigh- 
bor coupling possesses a refiection symmetry about k = 



7r/2. The next-nearest coupling does not have any re- 
flection symmetry except for fc = tt where it is odd so 
that no even 5 = two-triplon state can be excited. The 
spectral densities for both observables are dominated by 
the two-triplon bound state S'l which contains most of 
the spectral weight. This bound state exists for all mo- 
menta independent of dimerization and frustration. The 
binding energy increases by turning on the frustration. 

The spectral weight of the two-triplon continuum is 
concentrated at the lower band edge for all considered 
values of dimerization and frustration. For finite dimer- 
ization, the lower band edge displays a square root be- 
havior for a = 0.25 in accordance with the results of the 
sine-Gordon model. Divergences may occur only at ex- 
actly zero dimerization. The behavior changes similarly 
to the S — 1 case if a two-triplon bound state is almost 
degenerate with the lower band edge of the two-triplon 
continuum. Such a degeneracy appeared in our data for 
a ~ and for a ~ 0.5. 

Finally, we presented results for the Raman response 
and the infrared absorption. Both experiments are dom- 
inated for the values of dimerization considered by the 
bound state S'l. This bound state produces two van- 
Hove singularities in the infrared absorption resulting 
from k — 'k/2 and k — n. The van- Hove singularity at 
lower energies becomes more important for larger values 
of the frustration. 



VII. CONCLUSIONS 

We have shown that continuous unitary transforma- 
tions (CUTs) are an excellent tool to calculate spectral 
densities in high resolution for the dimerized and frus- 
trated Heisenberg chain. We provided a detailed study 
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FIG. 16: Two-triplon spectral density l2{k,uj) for -R^nn with 
A = 0.6 and a = 0.0 (a), a = 0.25 (b) and a = 0.5 (c). Gray 
lines denote lower and upper bound of the continuum. Black 
lines indicate dispersion of two-triplon bound states. 



of the spectral densities and extracted the generic fea- 
tures of the spectral properties. These data will help to 
analyze a large variety of spectroscopic measurements for 
dimerized and frustrated spin chains systems. 

We used a description in terms of triplons (elementary 
triplets) carrying S = 1. PreviousljiSS, we had shown 
for unfrustrated chains that a description in terms of two 
triplons is sufficient even in the limit of zero dimcrization. 
The triplon may also serve as an elementary excitation of 
the uniform Heisenberg chain besides the well-established 
spinon excitation. In the present work, we found further 
strong support for this result. The two-triplon spectral 
densities computed at finite dimcrization displayed well- 
developed precursors of the results for the uniform chain 
based on spinous, e.g., the dynamical structure factor 
probing the = 1 secto]4 or the optical absorption^^ 
probing the S = sector. 

In the frustrated case it is as yet an open issue 
whether a triplonic description works also in the limit of 
zero dimcrization. Especially the gapped phase (ao > 
Q^o,c) requires that the two-spinon continuum between 
A and 2A can be described by a dense distribution 
of bound many-triplon states. For this to occur an 
infinite-range effective interaction is necessary which is 
beyond the scope of the perturbative CUTs. But we 
take the nice agreement between the shifts of spectral 
weight obtained by complete exact diagonalization at fi- 
nite temperatures^ for the undimerized frustrated chain 
with our results at finite dimcrization as indication that 
the triplon description can be extended to the undimer- 
ized frustrated chain, too. 

The comparison of our results to those obtained by 
mapping the spin chain to a sine-Gordon model led to 
a number of insights. Both approaches agree that the 
generic singularity at the lower band edge is a square 
root, either a divergence or a zero. The divergence oc- 
curs if and only if a bound state is degenerate with the 
band edge. It turned out that the predictions of the sin- 
gle sine-Gordon model hold for critical frustration only 
in agreement with previous conclusions based on numer- 
ical results22,. The general spin chain at non-vanishing 
dimcrization requires to go beyond the single sine-Gordon 
model. We showed analytically that in the self-consistent 
harmonic approximation (equivalent to first order renor- 
malization) the Umklapp term is as large as the mass 
termS. 

The present study based on perturbative continuous 
unitary transformation is limited by the tractable max- 
imum range of hopping and interaction processes which 
correlates with the maximum order. Thus an investiga- 
tion of the spectral properties for small and zero gaps, 
i.e., dimerizations, will require to develop new methods 
allowing for larger ranges (correlation lengths). The con- 
tinuous unitary transformations do not need to be real- 
ized perturbatively. Ongoing research^ deals with a self- 
similar realization which renders the treatment of larger 
ranges possible. In this way, a closer look at critical sys- 
tems and systems with massive spinous will come within 
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FIG. 17: Two-triplon bound states for -Rnnn with A = 0.6 and a — 0.0 (a), a = 0.25 (b) and a = 0.5 (c). Left panels show the 
dispersion of the bound states; right panels the spectral weights of the bound states multiplied by the indicated factors. Gray 
lines denote lower bound of the continuum. 
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FIG. 18: Raman line shape for -Rnnn with additional broadening P — 0.01: (a) a = (b) a = 0.25 (c) a — 0.5. In each picture, 
curves for A = {0.3; 0.4; 0.5; 0.6} are shown, shifted for increasing A in j/-direction. The Raman response for a = is multiplied 
by 6. The insets in (b) and (c) zoom on the continua. 
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